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ABSTRACT 

The 85-foot radio telescope of the Naval Research 

Laboratory was used in an attempt to detect the 22,235 MHz 

transition (6_ 5 00 ) of H o 0 during the recent appearance 
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of Comet Bennett (1969i) . No H 2 0 emission of antenna 

temperature greater than 2.5 K was observed. We have derived 
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upper limits to the HgO column density (~ 2 x 10 molec/cm ) 
for various temperatures of the cometary gas. These limits 
have been compared with H 2 0 column densities calculated 
from two different cometary models. We find that on one 
of these models, our sensitivity was just at the threshold 
for detection of HgO. 
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Water is assumed to be the most abundant parent 
molecule in most proposed cometary models (Delsemme and 
Miller 1970) . There is no direct observational evidence 
for HO, but optical observations of the OH radical strongly 
imply the existence of a large quantity of HgO molecules 
near the comet nucleus. The recent detection of the 
6 - 5 00 rotational transition of H o 0 (Cheung et al. 

lO Zd 

1969) at a frequency of 22235.08 MHz (wavelength of 1.35 cm) 
in interstellar space provides a new impetus for searches 
for HgO in other astronomical objects. Comet Bennett (1969i) 
was a bright, young comet which reached its perihelion 
distance of 0.5 a.u. on 20 March 1970 at which time its 
geocentric distance was 0.7 a.u. An attempt was made to 
detect the 1.35-cm line or continuum emission from the head 
of this comet on 25 March, 28 March, 30 March, and 1 April 
1970. 

The 85 -foot (26-m) reflector at the Maryland Point 
Observatory of the Naval Research Laboratory was used. 

This antenna has a main-beam efficiency of 0.6 and a half- 
power beamwidth of 2.3 arcmin at a wavelength of 1.35 cm. 

An antenna temperature of 1 K for a point source is equivalent 
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to a flux density of 13 % 10 W m Hz . The beam- 
switched receiver employed a backward-wave oscillator which 
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was phase-locked to a tunable reference oscillator at 
350 MHz, which in turn was phase-locked to a quartz crystal 
standard for further frequency stability. The single-sideband 
system temperature was about 2500 K. The receiver gain was 
calibrated with an argon discharge tube of 200 K equivalent 
noise power and the receiver baseline was determined using 
"off-source" runs. Data were taken using both a bank of 10 
contiguous 100-kHz filters and a bank of 50 contiguous 
4-kHz filters. The 4-kHz profiles were averaged to an 
effective resolution of 12 kHz. A Doppler shift of 74.2 kHz 
is equivalent to a radial velocity of 1.0 km/sec. 

The primary limitation of this experiment was the fact 
that available ephemerides for the comet at the time of 
our observations were not of sufficient accuracy (within 
2 arcmin) for us to know the precise coordinates of the 
comet. The proximity of the comet to the sun precluded 
any attempts at alignment using an optical telescope. Thus 
we could only observe a grid of points centered on the 
"most likely" position of the comet and ex post facto 
determine which of those runs were taken while indeed ob- 
serving the comet head. In this way we were able to obtain 
useful data, but our integration times on the comet head 
were severely limited. 
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The results of our observations are summarized in 


Table I. Zone 1 is defined as a run in which the center 
of the beam was - 0.6 arcmin from the comet head and Zone 2 

i 

is a separation of 0.6-1. 2 arcmin. We can also place an 
upper limit of 10 K for the antenna temperature of any 
1.35-cm continuum radiation from the comet. 


TABLE I 

1.35-CM H 2 0 LINE SEARCH IN COMET BENNETT (1969i) 



Date 


Resolution 

(kHz) 

Zone 

Integration 
Time (min) 

Peak-- 

to- 

Peak 

ta(k) 

Velocity 

Range 

(km/sec) 

About 

V 

comet 

26 

March 

1970 

100 

1 

5 

2.5 

-8 -* +5 

1 

April 

1970 

12 

1 

20 

10 

-2 - +2 

1 

April 

1970 

12 

2 

15 

10 

-1 - +4 


I. 


The data in Table 1 may be used to estimate an upper 
limit for the column density of H^O in Comet Bennett. The 
observed antenna temperature T^ may be expressed as 

T A = Fri B T (1 " e ~ T) (1) 

where F is the beam dilution factor, r) the main-beam efficiency 

B 

(0.6), T the temperature of the cometary gas, and t the optical 
depth of the water vapor. Since the rotational levels 
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corresponding to the 1.35-cm line of II^O are 447 cm 
above the ground rotational level, the water vapor is 
optically thin for the cometary models considered in 
this paper (T - 500 K) . Equation 1 then reduces to 


T 


A 


7 'b t 


(t « 1) . 


where t is the optical depth of the water vapor averaged 
over the radio beam; for a rotational transition of an 
asymmetric molecule such as H 2 0, frequency-averaged over 
the width at half-amplitude, t is given by (Townes and 
Schawlow 1955) 


T 


8n 3/2 h 3/2 


3c k 


5/2 



-e/kT 

e 

^ 5 / 2 “ 


2 q 2 

p Sv 

Av“ 


<N> 


( 2 ) 


(3) 


where <N> is the mean column density of H^O in any state 

averaged over the beam response of the telescope, e the 

energy of the 6.,_ state, A, B, and C the rotational con- 

stants of the molecule, p the dipole moment, S the strength 

of the transition, v the frequency of the transition, and 

Av the half-width at half-intensity of the line. Substituting 

into Eqs. (2) and (3) the values A = 836 GHz, B = 435 GHz, 
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C = 278 GHz, p = 1.83 x 10 esu cm, S = 0.06, e = 0.0554 eV, 
v = 22235.08 MHz, and Av = (^£) ^ 2kT^tj ^“ , where m is the 
mass of the HgO molecule, we find 

, . -643/T 

T. = 1.2 x 10- 11 <N> - (4) 
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The minimum detectable antenna temperature was 2.5 K. 

The implied upper limits for mean HgO column densities 
are given in column 2 of Table II for comet gas tempera- 
tures ranging from 100 K to 500 K. 

Table II 

COMPARISON OF OBSERVATIONS AND MODEL CALCULATIONS FOR <N> 

, 2 

Mean Column Density <N> Molec/cm 


TOO 

Min. Detectable 
(T a = 2.5K) 

Fluid Dynamic 
Model 

Malaise Model 

100 

17 

13.0 x 10 

17 

0.6 x 10 

17 

8 x 10 

200 

2.1 

0.6 

8 

300 

1.6 

0.6 

8 

400 

1.7 

0.6 

8 

500 

1.9 

0.6 
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In order to interpret this lack of microwave emission, 

we compare our results with two cometary models: (i) the 

-2 

fluid dynamic model which predicts an r density distribution 
of molecules, and (ii) the model of Malaise (1970) derived 
from the observed behavior of the CN rotational distribution 
in several comets. 

In the fluid dynamic model, the density n„ as a function 

r 

of distance R from the nucleus is 

bER 2 

n F (R) = — t~ (5) 

vR 2 
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where E is the evaporation rate of the icy nucleus, R^ 

li 

the radius of the icy nucleus, v the flow velocity of 

i <\ 

s 

the gas, and b the number fraction of HgO molecules in 

the gas and is set equal to unity in the calculations 

which follow. It is assumed that HgO molecules only 

exist out to a radius R defined by the product of v and 

m 

the lifetime t of the H o 0 against photo-ionization. The 
m z 

expressions for integrating n (R) along the line-of-sight 

r 

and averaging over the radio beam have been derived by 

Huebner and Snyder (1970) for the fluid dynamic model. 

Column 3 of Table II gives the value of <N>„ using the 

£ 

18 2 6 
parameters b = 1, E = 10 molec/(cm sec), R = 5 x 10 cm, 

c 
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v = 10 cm/sec, and t = 1 x 10 sec (therefore R„ = 10 cm). 

m M 

The value for R c is an upper limit based upon the absence 
of gravitational perturbations due to comets on other 
bodies in the solar system. 

The empirical density distribution of Malaise, on the 
other hand, is given by 


V R) 


n R 
o o 


R 



( 6 ) 


where n Q is the gas density at the arbitrary distance R q 
q 

(taken as 10 cm) and R is the effective radius of the 

© 

molecule cloud, is a parameter fitted to observations which 
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varies between 5 x 10 and 2 x 10 cm, depending on the 

comet and its solar distance. The Malaise distribution 

does not strongly depend on the value of Il e , however, and 

g 

we adopt a value of 1 x 10 cm for our calculations. There 
are no reported observations of Comet Bennett that can be 

9 

used to determine n Q , but Malaise reports values of 5 x 10 
9 3 

and 3 x 10 molec/cm for Comet Seki-Lines (1962c) and Comet 

Ikeya (1963a) respectively. We adopt a value for n Q of 
9 -3 

4 x 10 cm for Comet Bennett since, although it was some 
2 I ?5 brighter than either Comet Seki-Lines or Comet Ikeya, 
its dust continuum was very strong relative to its gas 
emission. 

No simple exact analytic expression for <N>^, the 

integral of n (R) over the line-of-sight and the telescope 
M 

beam response, exists; however an examination of Eq. (6) 

shows that at worst a 30% error is introduced in the calcula- 

-R/R© 

tion of <N>., if (1-e ) is approximated as (R/R ) . Making 

this approximation, the column density N l.p) a projected 
distance p from the nucleus is 


N M (P) - 


2n R 2 cosh 1 (R /p) 
o o o 


R 


-R /R 
. o e 
1-e 


Finally <N>_, is obtained by integrating N_.(p) over the 
M M 


HgO cloud and averaging this value over the area of the 


M * 

(7) * 
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V 
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radio beam oi' radius R (R = 4 x 10 cm for these obser- 

B B 

vations) . The final expression is 


<n> m - 


2n R 
o o 


o -R /R 

R d 2 R (1-e ° e ) 

Be 


The calculated mean column density <N> 

2 17 
molec/cm as compared with 0.6 x 10 


model (Table II). 
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„ is 8 x 10 
M 

for the fluid dynamic 
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If all of the assumptions inherent is these calculations 
are correct, on the Malaise model our threshold of detection 
was almost adequate for a positive result. Column densities 
predicted on the fluid dynamic model, however, are an order 
of magnitude less. With only slightly improved sensitivity, 
a good radio test of the Malaise model may be possible for 
the next suitable comet. If these observations could be 
coupled with simultaneous optical observations, then the 
various model parameters would be more precisely known and 
more definitve statements concerning the mole fraction of 
H r O could be made. 

Ci 

The negative results of this and the formaldehyde search 
(Heubner and Snyder, 1970) both suggest however that under 
somewhat more favorable conditions i.e. smaller geocentric 
distance, larger telescope, higher sensitivity, longer observing 
times, the detection of parent molecules by radio observations 
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could be possible., assuming present ideas an nuclear 
composition are correct. A concentrated effort by radio 
astronomers on future nearly, bright comets is called for. 

We are indebted to B. Marsden who provided the latest 
orbits and ephemerides for Comet Bennett, L.W. Brown for 
assistance in the observing, and C.1I. Mayer for his support 
of this project. 
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